Integer quantum Hall effect and Hofstadter's butterfly spectra 
in three-dimensional metals in external periodic modulations 
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We propose that Hofstadter's butterfly accompanied by quantum Hall effect that is similar to 
those predicted to occur in 3D tight-binding systems by Koshino et al. [Phys. Rev. Lett. 86, 
1062 (2001)] can be realized in an entirely different system — 3D metals applied with weak external 
periodic modulations (e.g., acoustic waves). Namely, an effect of two periodic potentials interferes 
with Landau's quantization due to an applied magnetic field B, resulting generally in fractal energy 
gaps as a function of the tilting angle of B, for which the accompanying quantized Hall tensors 
are computed. The phenomenon arises from the fact that, while the present system has a different 
physical origin for the butterfly from the 3D tight-binding systems, the mathematical forms are 
remarkably equivalent. 



Introduction In two-diniensional(2D) periodic systems 
in magnetic fields, it is well-known that the interplay 
of Bragg's reflection and Landau's quantization gives 
rise to a fractal series of energy gaps, which is called 
Hofstadter's butterfly "'^ . Theoretically, the self-similar 
spectrum emerges both in the strong potential limit^ 
(tight-binding system) and in the weak potential limit^ 
(Landau-quantized system), where two cases relate to 
each other in their mathematical expressions. When the 
Fermi energy is within a gap in the butterfly the system 
exhibits the quantum Hall effect, where the Hall conduc- 
tivity is expressed as a topological quantum number for 
each energy gap^ . In recent experiments, a sigunature of 
Hofstadter's butterfly has been reported in the Hall con- 
ductivity measurement for semiconductor superlattices*. 

While the butterfly spectra had been known to be 
peculiar to 2D, we have previously found that three- 
dimensional(3D) tight-binding lattice can have a frac- 
tal energy spectra like Hofstadter's when the system is 
anisotropic (quasi-lD), where the 3D butterfly is shown 
to be a genuinly 3D effect rather than a remnant of a 2D 
butterfly^. This leads to an interesting transport phe- 
nomenon — quantum Hall effect (QHE) in three dimen- 
sions. According to the general argument^"^, QHE oc- 
curs even in 3D provided that there is an energy gap and 
that the Fermi energy lies in the gap, where the three 
components of the Hall conductivity tensor should be in- 
dividually quantized. 

Then an intriguing question arises: can we have a coun- 
terpart in the opposite limit, 3D weakly modulated sys- 
tems? Here we investigate the energy spectra in 3D con- 
tinuous systems with weak periodic potentials, and show 
that a modulation composed of two plane waves (such as 
two acoustic waves) gives rise to Hofstadter's butterfly, 
although the mechanism is quite distinct from one in the 
tight-binding system. We also find a clear mathemati- 
cal relationship connecting two limiting cases similarly 
to 2D, but we will see that three-dimensionality is essen- 
tial throughout. 

The problem has another interesting aspect, i.e., the 
QHE topological numbers. Here we have calculated the 



Hall conductivity in 3D weak potential systems, and 
find that the expression is remarkably equivalent to that 
for the corresponding energy gap in the tight-binding 
system^, which is unlike the relationship in 2D between 
the limits of weak ^ strong potentials. 

The present experimental setup (just an ordinary 3D 
metal with two acoustic waves applied) is so simple that 
this can be a strong candidate for experimentally detect- 
ing the 3D QHE and Hofstadter's butterfly. We shall 
estimate the required magnetic field and acoustic wave 
length, etc, for a possible realization in a semimetal sub- 
jected. 

Formulation: One modulation We first look at how a 
single modulation interferes with Landau's quantization 
due to a magnetic field B in a 3D uniform electron gas. 
The Hamiltonian is simply H = {l/2m){p-\-eA)'^ -i-U{r), 
where U (r) is a perturbative periodic potential, and the 
vector potential is taken as A = (0, Bx, 0) for B = 
{0,0, B). The eigenstates in the absence of U{r) are, as 
usual. 



\n, ky,kz) = Nn exp(z/cj^j/ -f ikzz)-^n{x/l + kyl). 



(1) 



where Hn is the Hermite polynomial with the Landau in- 
dex n, Nn a normalization factor, and I ~ y^fi/{eB) is the 
magnetic length 80 Afor _B = 10 T). The eigenenergy 
is En^k^ = {n-\- l/2)TiuJc + TT?'k1/{2m) with the electron 
mass m and the cyclotron frequency ujc = eB/m, which 
suggests that a 3D uniform electron gas becomes ID in 
that the motion along the magnetic field remains free. 
So, if we have an external modulation on top, we imme- 
diately expect energy gaps to emerge in the ID band. 
When we apply a single periodic potential U{r) = 
cos(G-r) in 3D, a mixing between the states \n, ky, kz) 
and \n' , ky ± MGy, k^ ± MG^) gives rise to energy gaps 
in the M-th order perturbation. The energy gaps due 
to the mixing between higher Landau bands will gener- 
ally tend to be hidden by other bands, so we focus on 
the lowest two levels. Namely, energy gaps should open 
within the lowest Landau level (n = n' = 0) and between 
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the lowest and second levels (n — O.n' = 1) (see, Fig. 
1). The corresponding matrix elements in the first order 
perturbation are 



(0, ky±Gy,k,±G,\U{r)\0,ky,k,) 
(1, ky±Gy,k,±G,\U{r)\0,ky,k,) 

= [;^e-«i'V4eT<fe.±^)G«;= ^ ^(G, - iGy)l/V2 (3) 



(2) 



where G± = yG^ + Gy. Since the matrix elements 

scale as exp(— the magnitude of the matrix el- 
ements, hence the size of the energy gap, becomes sig- 
nificant when Gl ^ 1 (i.e., exp(— G'"^/^/4) vanishes like 
exp(— const. /S) for B ^ 0). The condition Gl ^ 1 also 
suggests that the main energy gaps open for E < huJc, 
i.e., around the lowest or the second lowest Landau levels. 

In 3D systems, the Hall conductivity has three compo- 
nents. (T = (cTj/z, o"2;,.. CTf.y), where the Hall current in an 
electric field E is given by j = a x E. We can immedi- 
ately calculate cr for the present system when the Fermi 
energy Ep is in each of the energy gaps. If we consider 
generally the M-th order gap (M = 1,2,...) within the 
lowest Landau level (n = n' = 0), the number of states 
below Ep is 



Nf = 



MG^ 



eV 



2TrP 2Trh 



MGB, 



(4) 



where V is the system volume. From Widom-Stfeda's 
formula-''°'-^^, cr = ~{e/V){dNF / dB), we readily obtain 



^^ = -2^^^ (n = n' = 0). 
Similarly, we have for the inter-Landau gap 



2TTh 



2MG (n = 0,n' = l). 



(5) 



(6) 



A note is due here. 3D continuous systems having 
a one-dimensional modulation in magnetic fields have 
been investigated for a long time in terms of density- wave 
(DW) instabilities. Specifically, Halpcrin^ has calculated 
the Hall conductivity when the system is in a DW state, 
which corresponds to our calculation for the energy gap 
within the lowest Landau level. The inter- Landau level 
gap formation is new to the best of our knowledge. 

Two modulations Now, if we superpose two modula- 
tions Ui{r) = U\ cos(Gi • r) and U2{r) = U2 cos(G2 • r) 
having different wavevectors, the discussion above can 
be extended, where we can substitute MG with MGi + 
NG2 for the energy gap corresponding to Mth and A''th- 
order perturbations in J7i and U2, respectively. The Hall 
conductivity for the gap with the order M, N in the low- 
est Landau band thus becomes 



<T = - — {MGi+NG2) (n = 0,n' = 0), (7) 



where {M, N) becomes (2M, 2N) for the (M, iV)-th inter- 
Landau level gap. 

We can immediately realize that the spectrum should 
be sensitively affected by the commensurability of the 
ratio between the z-components, Giz and 6*22, since, for 
instance if the ratio is irrational MG\z + NG2Z can take 
continuous values and the spectrum should have gaps 
everywhere in the energy axis. If the magnetic field is 
rotated in the plane spanned by Gi and G2, the energy 
spectrum is then expected to have a fractal structure 
since the commensurability between Giz and G2z varies 
in a complicated manner. 

The spectrum plotted in Fig. 2 against the tilting an- 
gle of B obtained numerically, has indeed butterfly-like 
structure in the lowest Landau band as in Hofstadtcr's 
spectum. The Hall integers M,N have been calculated 
rigorously for each gap, by tracing the number of states 
below the gaps for the tilted field. In higher energy re- 
gion where two Landau levels overlap {E > {3/2)hcOc) 
we see the gap structure due to the inter-Landau level 
mixing (the gap labeled with M,N = 4,0(0,4) and 2,2 
in the figure), as well as the intra- level gaps which are 
visible only when gaps within n = and n = 1 happen 
to overlap. In the figure, the gap 3,1 (1,3) is a composite 
of 2,1 (1,2) in the lowest level and 1,0 (0,1) in the second. 

Correspondence with the strong potential case As men- 
tioned, we have previously found a butterfly-like spectra 
for the first time in periodic 3D systems which is modelled 
by the anisotropic (quasi-lD) tight-binding lattice. There 
we have shown that the period along the most conduc- 
tive direction (x) is not responsible for the emergence of 
the 3D butterfiy, while other two periods (along y, z) are 
relevant. Thus the two (present and the previous) cases 
are both 3D system with two periods and the only differ- 
ence is the strength of the periodic potential. One might 
then be tempted to think that the two butterfiy spectra 
cross over to each other when the amplitude of the peri- 
odic modulation is increased or decreased. However, as 
shown in the following, they have distinct physical ori- 
gins, residing in different limits, and the resemblance in 
the spectra comes from a mathematical relationship. 

To show this, let us concentrate on weak modulations 
with Gi,G2 ^ 1// for the lowest Landau level (i.e., 
low-energy region of the spectrum). To clarify the re- 
lationship, we concentrate on a case where B lies on the 
plane spanned by Gi, G2. If we take Gi and G2 on zx- 
plane (with Giy = G2y = 0) so that ky is conserved, 
Schrodinger's equation becomes 



2 7.2 



2m 



(8) 



where ipikz) is the amplitude of \0,ky,kz) with ky being 
constant, Uk,,k', is cq.(2) summed over Ui and U2 with 
Gy = 0, and is a function of ky. By Fourier-transforming 
with respect to z we obtain 
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■ijj{z)+Ui cos{GizZ - kyGixl'^)4>{z) 



+U2 cos{G2zZ - kyG2j^)i^{z) = E^iz), (9) 



(T = 



where Uj = C/jC"'^?^ with Gj_l = 
1,2). 

Now let us recapitulate how a butterfly spectrum arose 
in our original model on a 3D tight-binding lattice^. We 
consider the 3D orthorhombic lattice with the nearest- 
neighbor transfers t^, ty, along x, y, z, respectively. Wc 
assume here that the magnetic field B is applied parallel 
to y2;-plane {B^ = 0), which corresponds to the assump- 
tion made in the weak potential case that Gi, G2, B arc 
co-planar. Wc also assume that the system is quasi- ID 
(t^ ^ ty,tz), and apply the effective mass approxima- 
tion to the motion along the conductive direction x, to 
obtain the low-energy spectra (< ty or t^ from the bot- 
tom). In a gauge A = {0,BzX, —Byx) we can take the 
basis \E'(r) = e'^'^w+^'^'^xpix), and Schrodinger's equation 
becomes one-dimensional^'^: 



2m dx^ 



V'(x) 



'2ty COS 



eB,b 



kyb ] ij{x) 



—2tz COS 



pR r 

-^^x + k,c]ij{x)=Ei;{x), (10) 



where b, c are the lattice constants along y, z, respec- 
tively. Now we can see that eq.(9) and (10) have the idc;n- 
tifal form: ID equation with a double period. We should 
recall that Hofstadter's butterfly emerges generally in the 
spectrum of a doubly periodic ID system when the spec- 
trum is plotted against the ratio of the two periods^. 
In each of two equations, the ratio changes continuously 
when the magnetic field is rotated relative to the peri- 
odic potential, which is why Hofstadter's butterfly arises 
against the tilting angle in both cases. 

From the correspondence we can also note that t/i, C/2 
in eq.(9) correspond to ty,tz in (10), respectively, but 
their physical meaning is opposite in the following sense. 
Namely, Ui,U2 — » describes the limit of weak peri- 
odic potential, while ty,tz — > the strong limit. The 
equations become purely ID in both limits, where the 
system reduces to 'Landau tubes' (ID motion along B 
and the cyclotron motion on the plane _L B) in the for- 
mer, while in the latter wires along x confined by the 
strong potential on yz plane. On the other hand, if the 
perturbations {Ui,U2 or ty,tz) are too large, the mixing 
between different ID channels (i.e., different Landau lev- 
els in the former, different bound modes on yz plane in 
the latter) becomes so strong that the one-band approx- 
imation breaks down. We thus conclude that the two 
butterfly spectra in fact correspond to opposite limits 
(weak/strong potential). 

We can also establish a relationship between the Hall 
conductivities between the two cases. In the strong po- 
tential the Hall conductivity for the gap with M-th in 
ty and the A'^-th in tz in the tight-binding band can be 
obtained similarly as^ 



f N M 



(11) 



Since the primitive reciprocal lattice vectors are Gi = 
(27r/c)ez and G2 — (27r/6)ej^, this coincides with the 
corresponding expression in the weak potential, eq.(7). 
This is rather remarkable, since the corresponding wave 
functions in the two cases have totally different spatial 
behaviors, while they carry the identical Hall current. 
This contrasts with the 2D case, where the Hall integers 
are different between the corresponding gaps in the weak- 
and strong-potential limits. 
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FIG. 1. The energy dispersion plotted against the 
wavenumber along z (|| B) in a continuous 3D metal in 
a magnetic field B applied in a general direction (inset). 
Dashed arrows indicate where the modulation having Gz (the 
^-component) makes the dispersion gapped. 



Some notes on the quantum Hall effect are due here. 
First, the non-monotonic behavior in the Hall conduc- 
tivity as Ep is increased (see. Fig. 2) should be a hall- 
mark of the 3D butterfly spectrum as in the 2D butterfly. 
As for the quantum Hall plateaus, we have to consider 
the localization effect in the presence of disorder. While 
this is an interesting future problem, we speculate that 
each subband would evolve into the localized and the ex- 
tended states with mobility edges since we have a 3D 
system here, and that the Hall conductivity would be 
constant as long as Ep stays in the localized region. For 
the usual 2D butterfly, a numerical study for a dirty, fi- 
nite system-^^ shows that we still have a nonmonotonic 
behavior as a sign for the butterfly when the disorder is 
not too strong, so we expect a similar behavior in the 
disordered 3D butterfly as well. 

Second, while our calculation of the Hall conductivity 
is based on the bulk description, a finite system has edge 
states in the bulk gaps, which also contribute to the Hall 
current. We have previously shown^^ that the Hall con- 
ductance in a finite 3D system is still quantized, in the 3D 
QHE condition, when we take into accoimt both of the 
contributions from the surface current (which we have 
called the wrapping current) and from the bulk current, 
where the quantized values exactly coincides with those 
for the infinite system. The argument is quite general, 
and applies to the present problem as well. 
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FIG. 2. The energy spectra for 3D continuous systems in 
two periodic modulations having wavovcctors Gi. Gz (top), 
along with the plot for one modulation Gi (bottom), for 
\Gi\ = a/2//, Gi -L G2, and the amplitudes of the periodic 
potentials of 0.2'hciJc- The horizontal axis is the angle of B 
rotated in the plane spanned by Gi and G2. A pair of inte- 
gers (M, N) attached to each energy gap represents the Hall 
integers, eq.(7). The spectrum has been plotted for rational 
angles, tan^ = p/q {p,q: mutually prime integers), with an 
energy cut-off at 4.5fiu)c. The arrows on the right indicate 
the energy regions for the Landau bands in the un-modulated 
case. 



Experimental feasibility Finally let us comment on the 
experimental feasibility for the weakly modulated QHE 
system proposed here. To observe the QHE, Ep have 
to reside in an energy gap, which appears only in the 
lowest Landau levels. So we require the situation where 
only a few levels are occupied in 3D (the quantum limit), 
for which we need a large magnetic field or a small elec- 
tron concentration. We expect that semimetals should 
be suitable for the latter condition. We can estimate the 
required magnetic field for the quantum limit in bismuth, 
for example, to be B > lOT, which is quite modest. More 
stringent is the condition for the modulation, since the 
wave length should be such that the Fermi level is in the 
gap created by the modulation. If we take the acoustic 
wave for the periodic potential and plug in the velocity 
of sound for Bi, we have / ~ lOOGHz for B = lOT, which 
is rather high, although the required frequency (oc Fermi 
wavenumber) decreases with B like / oc 1/B. The cou- 
pling of the electrons with the acoustic wave, on the other 
hand, is dominated by the piezoelectric coefficient of the 
material, so compounds may be advantageous in this re- 
spect. An obvious advantage of the externally modulated 
system proposed here over the tight-binding system is 
that we can change the wavelength of the external mod- 
ulation at our disposal. So we expect that we can find 
a wider possibility for observing Hofstadter's butterfly 
with the 3D QHE. 
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